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摘 要 : 基于 长 时 间 序 列 遥 感 数据 ,提取 能 够 反映 生态 环境 质量 的 4 个 关键 指标 ,利用 主 成 分 分 析 方 法 构建 了 遥 
感 生态 指数 (Remote Sensing Ecological Index, RSED ,结合 Mann-Kendall 趋势 显著 性 检验 分 析 内 蒙古 呼伦贝尔 草原 
1990—2021 年 生态 环境 质量 变化 情况 及 变化 趋势 。 结 果 表 明 :(1) RSEI 能 够 全 面 反 应 内 蒙古 草原 生态 环境 质量 。 
30a 间 草原 生态 环境 质量 以 改善 为 主 , 整 体 呈 现 每 年 0.0037 的 上 升 趋势 (P<0.01) ,线性 回归 尼 达 0.39。 通 过 M-K 趋 
势 显著 性 检验 ,RSEI 在 2012 一 2021 年 间 呈 现 显 著 上 升 趋势 ,遥感 生态 环境 质量 得 到 明显 改善 (2) 总 体 而 言 ,研究 
区 内 RSEI 的 植被 绿 度 对 第 一 主 成 分 的 贡献 最 大 ,因此 二 者 结果 存在 较 高 相似 度 ; 但 生态 环境 质量 的 构成 包括 绿 度 、 
FE 湿度 和 热度 等 多 要 素 ,单一 绿 度 难 以 完全 表征 区 域 生 态 质量 情况 .(3) RSEI 呈 现 出 由 西向 东 逐 渐变 好 的 空间 
分 布 状 况 。 研 究 区 西部 牧 业 草原 地 区 生态 质量 状况 差 于 东部 森林 农耕 地 区 ,RSEI 空 间 呈 现 趋势 与 植被 覆盖 度 相 
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同 。 本 研究 丰富 了 内 蒙古 草原 生态 环境 质量 评估 方法 ,可 为 内 蒙古 草 向 草原 的 生态 保护 .生态 修复 .开发 利用 和 持 


续 发 展 提供 重要 参考 。 


关键 词 : 呼伦贝尔 草原 ; 陈 巴 尔 虎 旗 ; 遥感 生态 指数 ; Landsat; 时 空 分 析 


生态 环境 质量 状况 与 人 类 生活 息息相关 ,协调 
草原 生态 健康 与 生态 环境 保护 间 的 关系 ,是 推动 我 
国 草 原 可 持续 发 展 的 核心 问题 之 一 。 正 确认 识 和 
评价 草地 生态 环境 质量 对 提高 草地 资源 利用 效率 
和 生态 环境 保护 具有 重要 的 现实 意义 "1。 

呼伦贝尔 草原 地 处 蒙古 高 原 东 部 草原 腹地 , 拥 
有 纯 天 然 草 多 草原 。 作 为 我 国 典型 的 气候 变化 敏 
感 区 ,呼伦贝尔 草原 的 生态 环境 质量 近 几 十 年 来 一 
直 受 到 气候 变化 和 人 类 活动 的 影响 中。 研究 表明 ， 
Jr 30 a 来 ,呼伦贝尔 市 陈 巴 尔 虎 旗 的 平均 气温 呈 明 
显 上 升 趋势 ,而 降水 量 呈 明显 下 降 趋 势 ,导致 呼 伦 
贝尔 草原 自然 灾害 频 发 ,给 草原 地 区 经 济 与 社会 可 
持续 发 展 带 来 了 巨大 威胁 站 。 由 于 内 蒙古 草原 生态 
环境 的 时 空 变化 对 内 蒙古 草原 未 来 资源 保护 、 开 发 
和 利用 具有 重要 影响 ,因此 有 必要 结合 遥感 时 空 数 
据 对 草原 生态 进行 更 加 精确 的 研究 ,探求 草原 生态 


收 稿 日 期 : 2022-05-28; 修订 日 期 : 2022-11-24 


的 变化 趋势 和 特征 ,为 进一步 改进 和 强化 内 蒙古 草 
原 的 资源 保护 和 合理 开发 利用 提供 科学 依据 与 数 
UB 

遥感 生态 指数 (Remote Sensing Ecological In- 
dex, RSET) 是 综合 评价 生态 环境 质量 的 有 效 方法 ， 
是 专门 为 利用 遥感 技术 评价 生态 状况 而 开发 的 一 
种 综合 评价 指标 “。RSEI 是 多 因子 变化 综合 分 析 ， 
较 单 因子 分 析 而 言 ,多 因子 综合 分 析 更 加 全 面 和 准 
确 。 而 且 ,RSEI 可 以 将 反映 生态 问题 的 绿 度 指标 
(NDVI) 、 湿 度 指标 (Wet) . 干 度 指标 (NDBSIT) 和 热度 
指标 (LST) 通 过 主 成 分 分 析 方 法 进行 耦合 ,形成 一 
个 可 以 量化 的 生态 环境 质量 指数 。 近 些 年 ,RSEI 分 
析 主 要 侧重 于 长 时 间 序 列 如 20~30 a 间 生态 环境 质 
量 状况 ,空间 上 集中 于 如 矿区 、 城 市 .流域 .村 落 、 海 
岸 线 等 ””。 如 排 日 海 . 合 力 力 等 “以 新 疆 乌 鲁 木 齐 
市 为 研究 区 ,利用 RSET 对 当地 生态 环境 进行 了 评 
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价 ;高 鹏 文 等 "为 了 了 解 哈密 市 生态 环境 变化 ,对 干 
旱 区 城镇 哈密 市 伊 州 区 的 生态 质量 进行 遥感 动态 
监测 ,研究 结果 表明 生态 环境 质量 总 体 向 好 ; 张 亚 
球 等 "使 用 RSEI 对 区 县 级 生态 环境 进行 评价 研究 。 
GEE 利用 Google 的 海量 计算 能 力 可 以 解决 各 
种 具有 重大 影响 的 社会 问题 中 。Gorelick 等 在 
研究 中 将 GEE 数 据 .系统 架构 数据 分 布 模式 、 效 率 
性 能 .应 用 与 未 来 挑战 都 进行 了 详细 介绍 。Lu 等” 
基于 GFE 平 台 和 中 分 辩 率 成 像 光 谱 仪 (Moderate-res- 
olution imaging spectroradiometer, MODIS) 数 据 ,利用 
RSEI 分 析 了 北京 近 20 a 来 生态 环境 质量 的 时 空 变 
化 ,研究 表明 ,2001 一 2020 年 北京 生态 环境 质量 有 
所 改善 。 因 此 ,考虑 到 目前 内 蒙古 草原 生态 环境 研 
究 的 方法 和 手段 还 相对 较 少 ,对 其 生态 系统 和 环境 
质量 变化 的 研究 还 缺少 完善 的 体系 ,本 文 利用 与 生 
态 环境 质量 有 关 的 气候 条 件 、 大 气 环境 .植被 状况 
等 遥感 数据 ,选用 与 生态 环境 质量 直接 相关 的 4 个 
指标 ( 绿 、 湿 \ 干 、 热 ) ,结合 GEE 平 台 构 建 内 蒙古 草 
原生 态 环 境 时空 变 化 综合 指数 ,重点 对 内 蒙古 牧 业 
草原 放牧 + 农田 耕作 的 草原 生态 环境 变化 情况 进行 
系统 探究 ,对 内 蒙古 草原 生态 文明 建设 具有 重要 


1 数据 与 方法 


1.1 研究 区 概况 
研究 区 位 于 内 蒙古 自治 区 呼伦贝尔 市 陈 巴尔 
虎 着 ,属于 呼伦贝尔 草原 核心 区 域 ' 引 ,如 图 1 所 示 。 
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c3 呼伦贝尔 市 


40?N 45?N 50?N 55?N 


35?N 


100°E 110°E 120°E 


MAP 48*48' -50*12'N, 118?22'-121*02'E, WMX K H 
面积 约 1.86x10' km ,平均 气温 在 0 左右 ,属于 典型 
的 中 温带 大 陆 性 草原 气候 ,四 季 气 候 特 点 明显 “。 册 
东北 向 西南 地 势 逐 渐 降 低 , 人 研究 区 内 植被 类 型 主要 
为 森林 和 草包 草原 两 类 "”。 

1.2 数据 来 源 与 方法 

1.2.1 过 感 数据 来 源 使 用 的 卫星 还 感人 资料 由 美国 
地 质 调查 局 (United States Geological Survey, USGS) 
(https://earthexplorer.usgs.gov/) 提供 ,数据 源 为 
1990 一 2021 年 陆地 卫星 Landsat 5、Landsat 7 和 


Landsat 8 影像 ,空间 分 辨 率 均 为 30 m, 通 过 GEE 平 
台 进 行 数据 处 理 。 


鉴于 3 个 传感器 波谱 差异 ,考虑 到 数据 稳定 性 
与 可 对 比 性 ,其 中 Landsat 5 3H] 1990—2011 年 7 一 
9 月 的 反射 率 , 如 1990 年 选用 7 月 7 日 .7 月 23 日 、8 
月 17 日 .8 月 24 日 .8 月 26 日 和 9 月 2 日 六 期 影像 ; 
Landsat 7 选用 2012 年 7 一 9 月 的 反射 率 ; Landsat 8 
选用 2013 一 2021 年 7 一 9 月 的 反射 率 。 分 别 计算 每 
EI S B NDVI, WET, NDBSI fll LST 4 个 指标 ( 绿 、 
WEST GA), HB LSTARDHISCK BRL ER 3808 
标 为 中 值 合成 。Landsat 数 据 包括 主题 成 像 仪 TM、 
卫星 传感器 ETM+ 和 OLI3 种 数据 ,为 尽量 消除 卫星 
影像 获取 时 间 差 异 带 来 的 结果 的 不 确定 性 ,充分 体 
现 研究 区 植被 信息 ,遥感 影像 成 像 时 间 选 定 为 7 一 9 
月 草原 生长 最 旺盛 阶段 "% 。 数 据 异 常年 份 采 用 前 
后 2 a 数 据 平均 值 蔡 代 ,对 所 有 数据 均 进 行 了 数据 中 
值 平滑 处 理 。 
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注 : 底 图 采用 自然 资源 部 标 椎 地 图 制作 , 审 图 号 为 GS(2016)2556 号 ,对 底 图 边界 无 修改 。 下 同 。 
图 1 研究 区 概况 
Fig. 1 Study area 
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12.0 iÉ ASHA 
(1) 绿 度 指标 
绿 度 指标 用 归 一 化 植被 指数 (Normalized Differ- 
ence Vegetation Index, NDVD 来 代表 ,计算 公式 为 : 
NDVI =(P yr 7 Pu (Pos +Prea) (1) 
式 中 :pwn mm 分别 代表 遥感 影像 近 红 外 波段 (Land- 
satS 遥 感 数据 中 第 4.3 波 段 ) 和 红 波 段 (Landsat 8 中 
第 5.4 波 段 ) 的 反射 率 。 
(2) 湿度 指标 
缕 帽 变换 得 到 3 个 湿度 分 量 (Web , 即 亮度 RE 
和 湿度 。 其 中 湿度 分 量 反映 植被 .土壤 .水 体 中温 
度 值 ,是 生境 变化 在 生态 研究 中 重要 指示 指标 。 不 
同 卫星 传感器 的 表达 式 不 同 ， 
Wet, 2 0.0315p, + 0.2021p, * 0.3102p, + 
0.1594p, — 0.6806p, — 0.6109p， 
Wet, 20.1510, +0.1972p, +0.3283p,+ 1 
0.3407p, — 0.7117p, — 0.4559p, 
式 中 :Wet 表 示 卫 星 监测 湿度 指标 ,由 蓝 、 绿 、 红 、 近 
红外 和 短波 红外 的 地 表 反 射 率 计 算得 到 ;pi~ps 和 p， 
分 别 对 应 Landsat 5 的 第 1~5、7 波 段 和 Landsat 8 中 
第 2~7 波 段 。 
(3) 干 度 指标 
土壤 指数 (SD 和 建筑 指数 (BD 2 个 指数 合成 代 
表 干 度 指标 ,表达 式 为 : 
NDBSI=(SI+1BD/2 
SI 7 (Psw * Pn) 7 (Dna. + Pun (Dsum * Dn.) + 


(2) 


(4) 
(5) 


(Prue tu) 
iL. 2Ü enu Dun P oed 
IBI 
Ng tpa) | (Prin +Prea) d (Pe ed ts 
(6) 
2pswm Dun Pc ;reed 
十 十 

Los * p) ui *pya) (Po reg 1 Ts) 


式 中 :NDBSI 为 干 度 指标 ;IBI 为 建筑 指数 ;SI 为 土壤 
指数 ;pu posa Dis Dun svi 2] 31] 7J Landsat 5 中 第 
1、2、3、4、5 波 段 反 射 率 和 Landsat 8 中 第 2、3、4、5.、6 
波段 反射 率 (分 别 为 红 波段 . 绿 波 段 、 蓝 波段 、 近 红 
外 波段 和 短波 红外 1 波段 )。 

(4) 热度 指标 

选用 地 表 温 度 ,利用 卫星 监测 遥感 数据 热 红 外 
波段 辐射 值 加 最 新 修订 定 标 参数 来 代表 ,表达 式 为 : 

LST = T/[1 +(A T/p)lne] (7) 


其 中 : 


T=K,/n(K,/L,+ 1) (8) 
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式 中 :LST 为 热度 指标 ;7 为 传感器 温度 值 ;为 
Landsat 5 数据 的 波段 6, 中 心 波长 为 11.5 km, 或 
Landsat 8 数据 的 波段 10, 中 心 波 长 10.9 um; p X 
1.438 x 107mK;e 为 地 表 比 辐射 率 ,根据 Sobrino 提 
出 的 NDVI 阅 值 处 理 获得 ;K、K; 为 最 新 修订 定 标 参 
数 ;5 为 辐射 值 ;五 为 热 红 外 波段 的 增益 值 ;Zr 为 像 
元 灰 度 值 ; 环 为 偏 置 值 , 可 从 遥感 影像 头 文件 中 
获取 。 
1.2.3 技术 路 线 本 研究 的 具体 工作 流程 如 图 2 所 
示 。 首 先 , 基 于 GEE 云 平台 上 存储 的 Landsat 影像 
数据 ,计算 NDVI WET, NDBSI 和 LST 4 个 遥感 指 
标 ,并 利用 主 成 分 分 析 (PCA) 生 成 了 1990—2021 年 
陈 巴 尔 虎 旋 RSEI 的 时 空 分 布 图 。 在 此 基础 上 ,运用 
M-K 趋 势 分 析 方 法 对 陈 巴 尔 虎 旗 生 态 质量 变化 趋势 
进行 分 析 和 检验 。 
13 RSEI 遥 感 生态 指数 模型 构建 
采用 PCA 构建 RSEI, 对 每 个 主 成 分 的 贡献 度 主 
动 确定 ,将 NDVI.WET、NDBSI 和 LST 4 个 指标 信息 
集中 于 第 一 到 第 二 个 主 成 分 中 ,得 到 第 一 主 成 分 特 
征 向 量 PC1。 然 而 ,由 于 遥感 数据 提取 的 4 个 指标 
量 纲 不 同 ,每 个 指标 都 需要 归 一 化 处 理 。 同 时 ,为 
了 避免 大 量 水 域 面积 影响 主 成 分 分 析 载 荷 分 布 , 采 
用 水 体 指数 (NDWD 将 各 指标 中 水 体 掩 膜 处 理 后 , 进 
行 指标 归 一 化 和 主 成 分 分 析 。 指 标 归 一 化 方法 
如 下 : 


NI, =(1- 1,9, 7 Lu) (10) 
式 中 :MD 为 第 ;种 指标 归 一 化 后 的 值 , 取 值 范围 为 
[0, 1] ;7 为 各 个 指标 在 像 元 ;的 值 ; 心 ,为 该 指标 的 最 
小 值 ;为 该 指标 的 最 大 值 。 
RSEIL, = 1 - PCI(NDVL WET,LST,NDBST (11) 
RSEI = RSEI, -RSEIL /RSEL -RSEIL， (12) 
对 归 一 化 后 的 4 个 指标 进行 主 成 分 MM. 不 同 
年 份 分 别 进行 PCA 分 析 后 ,每 年 数据 提取 PC1 ,并 通 
过 1 减 去 PC1 获得 每 年 初始 的 RSEL, 以 保证 数值 越 
大 其 生态 质量 越 好 ;对 RSEl, 进 行 归 一 化 处 理 , 方 便 
生态 环境 质量 比较 。 最 后 ,运用 等 间距 分 级 方法 ， 
对 RSEI 每 0.2 划 分 一 个 等 级 , 共 5 个 等 级 ,分 别 代表 
生态 质量 差 较 差 一 般 、 良 好 和 优 。 
1.4 M-K 趋 势 显著 性 检验 
M-K 趋势 显著 性 检验 是 一 种 非 参 数 检验 ,用 于 
分 析 和 测量 非 线 性 系统 的 非 周 期 性 周期 长 度 。 不 需 
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图 2 技术 路 线 
Fig.2 Technology roadmap 


要 数据 呈正 态 分 布 ,对 时 间 序 列 中 的 异常 值 敏 感性 
较 低 ,可 从 某 个 方面 揭示 时 间 序 列 的 演化 规律 呈 。 
鉴于 M-K 趋 势 显著 性 检验 方法 可 以 确切 的 得 到 突 
变 点 ,因此 本 文 使 用 M-K 法 进行 生态 环境 质量 突变 
检测 。 

对 于 nn 个 样本 的 时 间 序 列 x, 构 建 一 个 序列 : 


i ls, 
sb "212, k=1,2,…,n) (13) 
zm 


0,x, Sx, 
式 中 :秩序 列 S, 是 第 ;时 刻 数值 大 于 7 时 刻 数值 个 
数 的 累计 数 。 


定义 统计 变量 : 
S, cS 
E P (14) 
|var[S,] 
其 中 : 
BEAG Li D vais- NE (15) 


E[S] 2j S, 的 均值 ,var[Si] 为 S, 的 方差 。 
给 定 显 著 性 水 平 a, 如 果 IUFil > Us, ZR Pl 

有 明显 的 趋势 。 参 考 这 个 方法 逆序 ,可 以 存在 UB: 
UB,=-UF,, (k=n,n—1,.…)]) (16) 
式 中 ;UF 为 一 条 变化 曲线 ;UB 表征 正确 的 逆序 序列 
的 发 展 趋势 。 如 果 UF 和 UB 曲线 的 交点 在 显著 水 
平 a 之 内 , 则 交点 的 具体 年 份 确定 ,说 明 该 年 份 的 参 
数 呈 现 突 增 状态 。 基 于 Python 实现 了 自动 制图 , 直 


2 结果 与 分 析 


2.1 RSEI 与 各 指标 统计 值 分 析 

通过 PCA 对 研究 区 进行 生态 环境 质量 分 析 ,如 
表 1 所 示 PC1 最 稳定 ,特征 值 贡献 率 远 大 于 其 他 3 个 
主 成 分 特征 值 贡 献 率 。4 个 年 份 中 PC1 的 特征 值 所 
占 比 例 均 在 58% 以 上 ,说 明 PCA 分析 法 效果 显著 ， 
对 通 感 生态 指数 计算 具有 代表 性 。 其 他 主 成 分 分 
量 值 忽 大 忽 小 ,说 明 其 对 生态 指标 的 信息 包含 不 完 
整 。PC1 中 NDVI 和 WET 为 正 值 ,说 明 NDVI 和 WET 
对 生态 环境 起 到 了 正面 效应 ,而 LST 和 NDBSI 指 数 
呈 负 值 ,说 明 其 对 环境 产生 了 负面 影响 ,符合 实际 
生态 情况 ,与 前 人 研究 一 致 ””” 。 

由 表 2 可 知 , 1990—2020 4E 38] [8] ,研究 区 内 
RSEI 呈 绥 慢 变 好 趋势 ,RSEI 均 值 由 0.46 上升 至 
0.52 ,总 体 上 升 了 约 13.13% ,30 a 间 RSEI 总 体 趋 势 
为 先 下 降 后 上 升 。WET 和 NDBSI 均 值 基本 保持 不 
变 , 但 2020 年 WET 均 值 低 至 0.283。LST 均 值 基 本 
呈现 逐渐 增 大 趋势 。NDVI 均 值 由 1990 年 逐渐 下 
降 ,至 2020 年 出 现 最 大 值 0.705。 总 体 说 明 虽 然 内 
蒙古 呼伦贝尔 草原 农 牧 结合 发 展 进程 快速 ,但 呼 伦 
贝尔 草原 坚持 生态 优先 .绿色 发 展 理念 后 ,人们 环 
境 保护 意识 的 不 断 提升 ,草原 植被 长 势 总 体 偏好 、 


观 展现 UF 和 UB 曲线 的 交汇 点 ,用 于 M-K 突 变 检 测 。 


草原 生产 能 力 明 显 提高 ” ,使 得 草原 畜牧 业 发 展 并 
没有 对 内 蒙古 呼伦贝尔 草原 造成 巨大 伤害 ””。 
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表 1 主 成 分 分 析 结 果 
Tab. 1 Principal component analysis 
1990 年 2000 年 
PCI PC2 PC3 PC4 PCI PC2 PC3 PC4 
NDVI 0.660 -0.601 -0.437 -0.113 0.64 -0.593 -0.488 0.0204 
LST -0.352 -0.775 0.524 0.0357 -0.611 -0.778 0.142 -0.046 
NDBSI -0.588 -0.135 -0.555 -0.573 -0.432 0.205 -0.827 -0.295 
WET 0.308 0.145 0.476 -0.811 0.177 -0.039 0.239 -0.954 
特征 值 0.0201 0.007 0.0028 0.0008 0.0335 0.0112 0.0036 0.0005 
特征 值 贡献 率 /% 65.50 22.85 9.12 2.59. 68.51 22.99 7.41 1.09 
参数 2010 年 2020 年 
PCI PC2 PC3 PC4 PCI PC2 PC3 PC4 
NDVI 0.764 0.517 0.384 -0.0401 0.549 -0.692 -0.468 0.0231 
LST -0.477 0.855 -0.203 0.00253 -0.331 -0.695 0.639 -0.015 
NDBSI -0.302 -0.0324 0.563 -0.703 -0.732 -0.190 -0.593 -0.277 
WET 0.313 0.00985 -0.768 -0.639 0.229 0.0488 0.15 -0.961 
特征 值 0.0309 0.0197 0.00239 0.000525 0.0331 0.0071 0.0035 0.0004 
特征 值 贡献 率 /9 57.77 36.78 4.47 0.98 75.05 16.05 7.95 0.95 
表 2 各 年 份 统计 结果 
Tab. 2 Statistical results of each year 
25 1990 4E. 20004E 
最 大 值 最 小 值 平均 值 标准 差 最 大 值 最 小 值 平均 值 标准 差 
RSEI 0.992 0.0186 0.463 0.164 0.992 0.0240 0.412 0.132 
NDVI 0.957 -0.597 0.610 0.109 0.973 -0.465 0.559 0.136 
LST 32.248 3.942 19.578 2.402 46.605 10.466 27.241 4.648 
NDBSI 1.600 -3.228 -0.0694 0.0978 0.632 -0.731 0.0414 0.104 
WET 0.923 -0.820 0.627 0.056 1.0044 -1.122 0.772 0.0419 
参数 2010 年 2020 年 
最 大 值 最 小 值 平均 值 标准 差 最 大 值 最 小 值 平均 值 标准 差 
RSEI 0.968 0.00632 0.481 0.179 0.994 0.00232 0.523 0.131 
NDVI 0.998 -0.696 0.552 0.160 1 -0.498 0.705 0.131 
LST 45.501 8.182 29.215 5.20 50.687 -19.225 37.574 6.424 
NDBSI 1.788 -11.760 0.0652 0.126 0.291 -0.729 -0.0371 0.132 
WET 0.99 -0.894 0.6005 0.0666 0.793 -0.227 0.283 0.0471 


2.2 土地 利用 类 型 变化 

研究 区 主要 土地 利用 类 型 为 农田 森林 、 草 地、 
水 体 等 几 类 3, 全 着 大 部 分 为 草地 ,东北 部 为 森林 ， 
农田 分 布 在 中 部 与 东部 区 域 。 图 3 中 草地 占 比 在 
81.7%~87.4% 左 右 , 和 森林 占 比 在 7.7%~11.3% 左 右 ， 
农田 占 比 在 3.8%~5.9% 左 右 ,水 体 占 比 在 0.52%~ 
0.69% 左 右 。1990 一 2020 年 土地 类 型 转移 中 草地 是 
转移 比例 最 高 的 土地 类 型 ,图 4 中 主要 列举 了 
1990 一 2020 年 间 转 移 最 多 的 3 种 土地 利用 类 型 , 草 
地 与 农田 类 型 的 转换 多 发 生 在 研究 区 的 东北 部 区 


域 。 从 表 3 中 可 以 看 出 ,1990 一 2020 年 间 草 地 覆盖 
与 其 他 各 类 土地 类 型 均 有 转 和 人 转 出 , 且 规 模 相 对 较 
大 ,土地 利用 结构 属于 不 稳定 状态 。 除 相同 土地 类 
型 转移 外 ,其 他 类 型 转移 中 农田 -草地 面积 约 
177.97 kn? , 占 农田 转 出 面积 的 26.89%; 草 地 -农田 
面积 约 550.47 km? , 占 草 地 转 出 面积 的 3.61%; 草 地 - 
森林 面积 约 590.34 km2 , 占 草地 转 出 面积 的 3.88%。 
2.3 RSEI 时 空 序列 变化 分 析 

图 5 为 陈 巴尔 虎 旗 各 年 份 RSEI 统 计 分 析 。 
2006 年 RSEI 平 均值 最 低 ,为 0.33, mi 2018 年 RSEI 
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图 3 2020 年 土地 利用 类 型 现状 
Fig.3 Land use type of 2020 
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图 4 1990 一 2020 年 土地 利用 类 型 转移 
Fig. 4 Land use type transfer from 1990 to 2020 


平均 值 达 30 a 最 大 值 (0.57), 较 1990 年 RSEI 平 均值 
(0.46) 增 长 了 0.11 ,涨幅 达 23.9%。NDVI 整 体 趋 势 
与 RSEI 走 势 表 现 相 近 ,其 中 ,2007 年 NDVI 平 均值 
最 小 为 0.49 ,最 大 值 出 现在 2021 年 为 0.76, 较 1990 
年 NDVI(0.61) 增 长 了 24.6%。LST 呈 现 增 加 -减少 - 
增加 的 波动 变化 ,波动 范围 在 20~40 % 之 间 。 因 研 
究 区 内 建筑 . 裸 土 分 布 相 对 较 少 ,所 以 NDBSI 平 均 
值 并 不 是 很 高 ,在 0.1 左 右 。 总 体 来 看 ,RSEI 均 值 在 


1990—2021 年 整体 呈 上 升 趋势, 表明 30 a 干 度 和 热 
度 对 生态 环境 的 负 效 应 明显 小 于 绿 度 和 湿度 对 生 
态 环 境 的 正 效应 。 

由 图 6 可 知 ,空间 尺度 上 ,由 于 研究 区 土地 利用 
类 型 不 同 ,研究 区 西部 分 布 多 为 草 多 草原 ,生产 方 
式 以 纯 牧 业 为 主 ;东部 多 为 农田 和 和 森林 结合 地 带 ， 
生产 方式 以 牧 业 + 耕地 为 主 ,植被 覆盖 度 呈 现 由 西 
向 东 和 逐渐 递增 的 趋势 "““。RSET 与 植被 覆盖 度 呈 现 
由 高 向 低 的 转化 趋势 一 致 ,可 以 清晰 地 观察 到 RSFI 
空间 分 布 格局 呈现 为 由 西向 东 逐 渐变 好 的 态势 ”。 
时 间 尺 度 上 ,从 1990—2021 年 ,草地 与 森林 、 农 田 
RSEI 值 在 120°E 附 近 出 现 明 显 界限 ,并 随时 间 增 加 
逐渐 东 移 。 
2.4 RSEI 等 级 变化 

由 图 7 可 知 , 各 年 份 RSEI 面 积 变化 除 优等 级 呈 
下 降 趋势 外 ,其 他 各 等 级 均 有 不 同 程度 的 上 升 趋 
势 , 其 中 较 差 等 级 趋势 变化 最 大 ,R=0.41, 其 次 为 一 
般 和 优 , 尼 分 别 为 0.35 和 0.34。 对 1990—2021 年 各 
生态 等 级 面积 进行 统计 分 析 , 优 占 比 13%~28.8% 之 
间 ,良好 占 比 25.1%~49.2% 之 间 ,一 般 占 比 20.3%~ 
36.896 , 较 差 占 比 1.5%~27.9% , 差 占 比 0%~2.3%， 
此 优 、 良 好 和 一 般 的 面积 占 比较 大 。 

由 图 8 可 知 , 陈 巴 尔 虎 旗 西 北方 向 草原 牧 业 区 
域 1990—2020 年间 RSEI 变 化 趋势 呈现 显著 变 好 ， 
其 他 大 部 分 是 略微 变 好 和 无 变化 部 分 。1990 一 
2000 年 , 2010—2020 年 间 大 部 分 为 略微 变 好 ， 
2000—2010 年 间 大 部 分 为 无 变化 与 略微 变 好 ， 
1990 一 2020 年 绝 大 部 分 是 略微 变 好 与 显著 变 好 , 表 
明 陈 巴尔 虎 旗 30 a 间 生态 环境 质量 状况 总 体 呈 现 变 


表 3 研究 区 1990 一 2020 年 土地 类 型 转移 


Tab.3 Land use type transfer in the study area from 1990 to 2020 /km 
1990 年 
草地 建筑 裸 地 农田 和 森林 湿地 水 体 总 计 
2020 Œ 草地 13993.36 0.057 11.00231 177.971 15.939 24.170 19.0455 14241.55 
建筑 29.073 2.17694 0.213 2.649 0.367 0.254 4.0404 38.774 
裸 地 48.155 0.00293 10.215 0.545 0 0.0351 0.902 59.854 
农田 550.467 0 0.0164 440.137 34.674 2.950 2.358 1030.604 
森林 590.342 0.0105 0.0240 36.940 1301.938 6.00475 15.216 1950.475 
湿地 0.114 0 0 0.0998 0.00704 1.188 0.0382 1.446 
水 体 17.476 0.109 0.127 3.598 0.0892 1.921 67.497 90.817 
总 计 15228.99 2.357 21.598 661.939 1353.014 36.522 109.096 17413.52 
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图 5 各 指标 时 间 序 列 变化 
Fig.5 Time series changes of each indicator 
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图 6 RSEI 等 级 
Fig.6 Ranking of RSEI 
好 趋势 ,其 中 2010 一 2020 年 间 变 好 趋势 最 高 。 平 ,在 此 期 间 RSEI 和 急剧 上 升 。 
2.5 M-K 趋 势 分 析 30 aH] , UF 曲线 与 VB 曲线 产生 1 个 交点 , 即 


如 图 9 可 知 ,从 1990 一 2001 年 间 RSEI 呈 现下 降 2010 一 2011 年 间 出 现 突变 点 ,但 范围 并 没有 突 
趋势 ; 自 2001 年 以 来 ,平均 RSEI 呈 现 出 较 明 显 上 升 ” 破 +0.05 置信 水 平 区 间 ,因此 可 以 初步 判定 陈 巴 尔 
趋 热 。UF HÆ 2012—2021 年间 远 远 超过 置信 水 “” 虎 旗 1990 一 2021 年 的 RSEI 变 化 整体 比较 平稳 。 
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14000[ .( 。 良 好。 一 般 7713781» 266101.— a 草 原生 态 环境 质量 得 到 明显 提升 。 该 研究 结果 与 
a 土地 生态 状态 类 似 ,植被 覆盖 度 与 土地 利用 类 型 对 

.. 10000 p rM cu Requests 3 研究 区 内 生态 环境 质量 都 起 到 一 定 作 用 2 。 

| 六 6 人 3028564 V196 86x-300613 由 图 10 可 知 ,1979 一 2021 年 间 呼 伦 贝 尔 草原 呈 

i a o 现 出 气温 升 高 的 发 展 趋势 ,使 呼伦贝尔 草原 阶段 性 

tapete 037l58961 气象 干旱 频 发 下。 对 气温 \ 降 水 与 RSEI 进 行 时 间 序 
2000 Sn 11578-23038 R?=0.34 、 i 
M eseh dae. - * s 列 分 析 发 现 ,1991 年 .2006 年 RSEI 较 低 ,但 2006 年 
SPN PD SDSS 于 a 平 均 E ; 2017 年 2018 年 RSEI a 
m 高 ,但 2 个 年 份 的 降水 量 均 低 于 30 a 平 均值 。 可 以 
看 出 ,气温 和 降水 在 一 定 程度 上 对 于 RSEI 起 到 一 定 
图 7 RSEI 等 级 面积 变化 的 综合 影 响 
Fig.7 Changes in the area of RSEI grades 草场 利用 方式 随 章 场 产 权 制 度 改 变 而 改 恋 
3 讨论 1983 年 ,呼伦贝尔 草原 开始 实施 经 济 体制 改革 , 推 


利用 RSEI 指 数 综合 评估 生态 环境 质量 状况 ,所 
得 结论 与 植被 覆盖 度 一 致 “。 植 被 覆盖 度 高 的 地 
方 生 态 环境 质量 较 好 ,植被 覆盖 度 低 的 地 方 生态 环 
境 质 量 相 对 较 差 ,说明 生态 环境 质量 综合 评估 了 植 
被 生长 状况 ,适用 于 对 植被 进行 整体 分 析 ,尤其 是 
对 不 同 植被 覆盖 度 和 土地 利用 类 型 。 与 前 人 研究 
一 致 ,可 在 一 定 程 度 上 反映 草原 畜牧 承载 能 力 、 反 
映 草 原 退 化 状况 方面 发 挥 重要 作用 1。 

通过 Mann-Kendall 突变 检验 分 析 , 陈 巴尔 虎 族 
RSEI 在 1990—2001 年 间 呈 下 降 趋 势 , 从 2001 年 开 
始 逐 渐 升 高 ,其 中 2012 一 2021 年间 显著 增加 , 近 10 
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行 草场 承包 到 户 政策 ,使 用 权 30 a 不 变 , 即 从 “游牧 ” 
到 “ 定 牧 ?利用 方式 的 演变 二 ,草场 利用 方式 在 一 定 
程度 上 影响 到 呼伦贝尔 草原 生态 环境 质量 。 为 保 
护 和 改善 草原 牧区 生态 环境 ,保障 我 国生 态 屏 障 安 
全 ,我 国 颁 布 实 施 了 众多 生态 工程 ,其 中 围栏 禁 牧 
是 一 种 可 以 增加 草地 生物 多 样 性 和 恢复 生态 服务 
功能 的 有 效 措施 号。 牧民 对 草地 资源 的 珍惜 ,积极 
履行 围栏 禁 牧 政策 ,使 得 1990 一 2021 年 间 呼 伦 贝 尔 
草原 植被 生态 环境 质量 呈现 出 了 先 短暂 下 降 、 再 度 
上 升 的 模式 ,总 体 呈 现 向 好 趋势 。 对 于 呼伦贝尔 草 
原 草 地 退化 的 影响 因素 深入 分 析 发 现 , 社 会 经 济 发 
展 因 子 对 呼伦贝尔 草原 退化 的 影响 更 大 ,其 中 陈 巴 
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图 8 RSEI 动 态 变化 监测 
Fig. 8 Monitoring diagram of RSEI level change 
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图 9 1990 一 2021 年 M- 多 趋势 分 析 
Fig.9 M-K trend analysis from 1990 to 2021 
一 气温 。 降水 (1) RSEI 可 以 对 植被 进行 整体 综合 分 析 ,尤其 
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图 10 气温 、 降 水 与 RSEI 时 间 序 列 分 析 
Fig. 10. Temperature, precipitation and RSEI time series 


analysis 


尔 虎 旗 草 地 退化 主要 影响 因素 是 总 人 口 数 ”。 草 
场 利 用 方式 改变 过 程 中 ,当地 牧民 对 保护 生态 环境 
关注 度 逐 渐 增 高 ,可 见 传统 牧 业 生 产 对 草原 植被 有 
一 定 影响 ,但 如 何 做 好 草 畜 平衡 状态 还 需 进 一 步 深 
入 人 研究 ,立足 于 生态 环境 可 持续 发 展 的 角度 上 ,资源 
循环 性 畜牧 业经 营 模 式 将 是 今后 发 展 的 重点 ”“。 男 
外 ,农业 生产 也 改变 了 草原 生态 植被 状况 ,因此 宜 
农 宜 牧区 应 很 好 地 界定 并 做 好 过 渡 带 的 修复 和 
保护 。 


4 结论 


气候 变化 和 人 类 活动 复杂 作用 影响 下 ,内 蒙古 
呼伦贝尔 草原 生态 环境 质量 正 发 生 改变 ,因此 利用 
多 源 遥 感 卫星 数据 分 析 其 1990 一 2021 年 生长 季 
(7 一 9 月 ) 生 态 环境 质量 ,得 到 以 下 结论 : 


是 对 不 同 植被 覆盖 度 和 土地 利用 类 型 ,1990 一 2021 
年 间 呼 伦 贝 尔 草原 植被 生态 环境 质量 呈现 出 先 短 
暂 下 降 再度 上 升 的 模式 ,总体 呈现 向 好 趋势 。 

(2) 人 研究 区 内 RSEI 的 植被 绿 度 对 第 一 主 成 分 
的 贡献 最 大 ,因此 二 者 结果 存在 较 高 相似 度 ; 但 生 
态 环境 质量 的 构成 包括 绿 度 . 干 度 .湿度 和 热度 等 
多 要 素 ,单一 绿 度 难 以 完全 表征 区 域 生 态 环境 质 


三 | 


里 o 


(3) 内 蒙古 自治 区 呼伦贝尔 市 陈 巴 尔 虎 旗 牧 业 
草原 森林 农耕 发 展区 域 生 态 环境 质量 在 30 a 间 绝 
大 部 分 优 于 放牧 区 域 ,RSEI 呈 现 出 由 西向 东 逐 渐变 
好 的 空间 分 布 状况 。 
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Detecting the change and trend of remote sensing ecological quality in Inner 
Mongolia from 1990-2021: A case study of Chenbarhu 
Banner of Hulunbuir City 


ZHANG Yusi", BAO Yuhai", HE Zhonghua’ 
(1. School of Geographical Sciences, Inner Mongolia Normal University, Hohhot 010022, Inner Mongolia, China; 


2. Key Laboratory of Remote Sensing and Geographic Information System of Inner Mongolia Autonomous Region, 
Hohhot 010022, Inner Mongolia, China; 3. Zhejiang Climate Center, Hangzhou 310056, Zhejiang, China) 


Abstract: The ecological civilization construction in Inner Mongolia, an important ecological barrier in northern 
China, is important for national ecological security, the stability of the border areas, and the sustainable economic 
and societal development of the minority areas. In this study, we extracted four key indices based on long-term 
satellite observations to describe the quality of the ecosystem. We constructed the remote sensing ecological 
index (RSED through principal component analysis. Furthermore, we adopted the Mann-Kendall (M-K) trend test 
for analyzing the changes in ecological quality in the Hulunbuir grassland of Inner Mongolia from 1990 to 2021. 
The RSEI can comprehensively reflect the quality of the ecological environment, and results revealed that the 
quality of grassland in Inner Mongolia has improved in the past 30 years, with an increasing trend of 0.0037 
every year (P < 0.01; R? = 0.39). The significant increasing trend of RSEI and the improved ecological quality 
from 2012 to 2021 was further verified through the M-K trend significance test. A single factor cannot fully 
characterize regional ecological quality so, four factors i.e., greenness, dryness, humidity, and heat were included 
in RSEI. The space quality showed a gradual improvement from west to east, with lower RSEI for grassland in 
the west compared with that for forest and farming areas in the east. A similar trend for quality was observed with 
the variation of the vegetation coverage. This study enriched the assessment methods of grassland ecological 
quality in Inner Mongolia, which can provide an important reference for the ecological protection, restoration, 
utilization, and sustainable development of meadow steppe in Inner Mongolia. 


Keywords: Hulunbuir grassland; Chenbarhu Banner; RSEI; Landsat; spatiotemporal analysis 


